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SUMMARY 



The dynamic. effects of alternative manpower policies and 
programs can be "pretested" in a cojnputer based simulation. It is 
generally recognized that the production of Doctorates depends to 
a large degree on the Doctorate-holding faculty . Because the 
doctorate holders are in great demand by the other sectors of the 
economy, a circular or a "feedback" situation exists. The problem 
is further cbmplicate'd by the availability of developed student 
talent and by various socioecbnomic conditions existing at different 
periods of time within two or three decades to the time a study is 
made. 



This Stiidy developed a conceptual and a mathematical model to 
study the production of Doctorate's, Masters and Baccalaureate 
degrees tod their feedtback into 'higher education. The model 
cottsisting of over 200 hoh-linear difference equations was programmed 
for cc^Uter 'simulatibh and validation against historical data. 
Simulation attempts to describe what haS happened in the past have 
thus far proved ih'cbh'cliisive i However, these failures are 
attributable to software problems encountered with the University of 
Wisconsin Burroughs B5500 computer installation. Currently efforts 
are under way to resolve these problems and the study is continuing 
albeit unsupported. Once the above phase is accomplished, the model 
can be used to prescribe what will happen in the future with a fair 
level of confidence. 



1 



o 



INTRODUCTION 



The output of higher education in general* and of the 
professional schools in particular* has been a subject of concern 
to various agencies of government* of foundations* both public 
and private* and to the professional associations involved for over 
half of a century. Data regarding the number of students enrolled 
and degrees granted have been diligently .collected, compiled, 
classified and reported. On the other hand, data regarding the. 
utilization of the graduates in academic positions, industry, 
government, etc., are much more sparse, although available. Data 
regarding the retention within our economy of foreign nationals 
who have successfully completed their schooling here is more sparse 
yet, or more precisely unavailable, even though foreign nationals 
comprise a large segment of the student body: in our professional 
and/or postgraduate programs of study. It is in this setting that 
planners are attempting to set national policy .concerning higher 
and/or professional education in order to satisfy anticipated or 
perceived needs for trained talent in the years ahead. 



Objectives of this. Study 



The objectives of thia study are: 

(1) To develop a simulation vehicle to pretest 
alternative manpower policies and programs . 

(2) To study the dynamics of the higher education 
sector within our society. 

(3) To study the dynamic interrelationship of the 
various subsectors of higher education in the 
rest of the economy. 

(4) To study the usefulness of currently available 
data and to indicate what additional data would 
be useful in determining long-range policy and 
planning . 

(5) To determine the impact on our economy in general 
(and the higher education sector in particular) , 
of the inflow and outflow of people from other 
nations . 
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(6) To study the possible impact of development in 
the technology of education upon manpower 
redistributions in the economy and the timing of 
same (i.e., increases in student-faculty ratios).. 



During the last two decades attention has been turned to the 
production and retention, within the higher educational sector of 
the economy, of people holding the doctorate degree. It . s 
generally recognized that the production of doctorates depends to 
a large degree on the doctorate-holding faculty. Because the 
doctorate holders, especially in the sciences, are in peat demand 
by the other sectors of the economy, a circular or a feed-back 
situation exists. The problem is further complicated by the 
availability of developed student talent and by various socio- 
economic conditions existing at different periods of time within 
two or three decades prior to the time a study is made. Various 
studies of these problems have been initiated by agencies of 
government and by privately financed foundations. 



In order to attempt to rectify an existing situation or 
provide for future needs, initiators of programs are in need of a 
rational methodology to Evaluate the effects of their program. 
Historical data and judgment have. not always yielded satisfying 
bases for decisions. Both the military and indust^ are recognizing, 
to an increasing degree, the need for quick and economical ways of 
evaluating the various effects of their decisions. For this they 
have turned to analytical methods. Computer-based simulation is 
one method successfully used to study the dynamics of complex and 
non-linear problems. 

Bolt, Koltun and Levine (1) made the first attempt, in 
to model and to study the dynamics of the problem at hand. Their 
work was Indeed landmark in nature, for it showed the advantaps 
of dynamic over statistical analysis of manpower probl^ while 
recognizing that the two are not mutually exclusive. The above 
model conceptually split the economy into two sectors , oi^ 
containing the entire educational establishment (Hi^er Education) 
and the other the rest of the economy. The flows of graduates 
between these two sectors were established recognizing the fact 
that Ph.D. holders do return to the academic sector to take up 
faculty or administrative positions while others leave for industry, 

government, etc. All of the flows including the so-called 

"feedback" flows were described by two linear difference equations _ 
which were simultaneously solved at various values of the equations 
parameters . 

This highly aggregated approach was considered unjustifiable 
by this investigator, in view of the inherent complexity of the 
problem studied and, in light of modem computer capability to solve 
large sets of differential or difference equations. Consequently, 
a less aggregated model was proposed (2). 
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This model broke the educational sector up into four segments; 
undergraduate programs, master’s programs, doctoral programs, and 
post-doctoral programs. It broke the other sectors of society 
employing college or university-trained people into segments 
according to the highest degree earned by those within the segment. 
Furthermore, it showed the retirement and other attrition sectors 
more or less as a sink outside of either of the above two sectors. 

The model delineated the various possible paths for population 
shifts between the segments. 

In a manner similar to (1) , (3) and (4) and to many works in 
the physical sciences, it accounted for all the net flows to a 
segment and the rate of accumulation of people within the segment. 

The equations recognized the fact that degrees, especially 
at the doctorate level, are not earned at a given time of the year 
throughout society. Some schools operate on the semester basis, 
others operate on a trimester, and yet others on the quarter system. 
When aggregated, an assumption of continuity in flows seemed a 
little more realistic than an assumption of discreteness.. Thus, 
the equations offered were differential equations. Depending on 
the postulates regarding lead- lag relationships in the production 
of the various degrees these equations were allowed t.p be 
non-linear. 

The basic advantages of that model over (1) were described in 
outline as follows: 

1) It recognized the input of students into the higher 
education sector. 

2) In the educational sector, it distinguished between 
persons who have recently become engaged in the educational 
function and persons who have worked in education for 

many years . 

3) It recognized that the relationships between the variables 
involved were not necessarily linear in nature. 

4) It distinguished between the use of doctorate holders in 
education at the doctorate, master's and undergraduate 
levels . 

5) It considered the effect of the rates of production of 
high school, bachelor and master's graduates in successively 
preceding years. 

6) It considered post-doctoral university programs and the 
Interrelated flows from them to teaching at the various 
levels and to the other sectors employing doctorates . 
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7) In as much as the r^te of doctorate production has changed 
drastically during the last three decades, because of the 
depression of the 1930’s, World War II, and the post-rvar 
and the cold-war periods, it emphasized that the age mix 
is not a uniform one. Thus, the number of degrees granted 
some years prior (say 30 or 35 years) is the independent 
variable used in calculating attrition in. that model. 

8) .Concern for economic, social, and physical influences in 

doctorate production and shifts in employment was. 
accomodated. 

A tangible result of this study was the rearrangement j 
extension and modification of the above model. This latter version 
of the model described in the fpllowlng section and spelled out, in 
detail in the Appendix to this report was discussed and reviewed 
with colleagues in various settings. The first presentation took 
place at the 14th international Meeting of the Institute of 
Management Sciences which was held in Mexico City, August 22-25 j 
1967. (6) The next presentation took place at the "Operations 
Analysis in Education" Symposium sponsored by the Office of Education 
and held in Washington, D.C., November 19-22, 1967. The Fall “issue 
of the journal, Soclo Economic Planning Sciences , whiCh is devoted 
to the Proceedings of thctt conference, will carry the paper ''On the 
Generatioxi of Doctorates and Their Feedback into Higher Education . 

A presentation of this model was also made at the "Engineering 
Manpower Planning Institute", which was held in Madison, Wisconsin, 
January 8, 9, 1968. At all of these meetings this model generated 
quite a bit of active discussion resulting in several additional 
changes of the model. | 

Unlike its predecessor, this model separates out the student 
and the faculty siibsectors within the academic sector and does this 
at all levels of educational attainment. It provides for inflows as 
well as outflows of foreign nationals at various levels of attain- 
ment.^ This model explicitly pirovides for the psychological, 
sociological and economic factors which influence the movement of 
personnel from one sector to another within our economy. It 
furthermore is expr'^.ssed in tsrms of ^^ifference" types of equations 
as opposed to differential equations. 



^This extension was suggested by Dr. Lindsey Harmon, Director 
of Research, Office of Scientific Personnel, National Academy of 
Sciences National Research Council. 

^The change from differential to difference equations was 
necessitated by the availability of the DYNAMO compiler for the 
B 5500 computer. Although difference equations are used, the time 
Increment considered is relatively short, specifically it is 1 
calendar .month . 
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A detailed description of this model is given in Appendix I 
to this report. 

RESULTS 

Because of the software problems with the Burroughs B 5500 
computer which was the only computer, large enough for the purpose 
available to us at the University of Wisconsin, this investigator 
regrets to report that no useful simulation runs are as yet 
available. Although the project grant has run its course,’ work is 
still continuing on the project, albeit at the principal 
investigator’s new location, the Department of Operations Research, 
Case Western Reserve University, Cleveland, Ohio. Two graduate 
students are working with the principal investigator on this 
project. Current work involves the simultaneous programming of 
this model for solution on either th4 Univac 1108 computer or the 
IBM 360-40 machines. It is the investigator’s intent to submit a 
supplementary final report on this work a& soon as the simulation 
has been completed resulting in substantive or new insights to the 
problem. 



METHOD; MODEL DEVELOPMENT 

Figure 1 describes all of the sectors and sub sept or s of our 
economy and the paths through which population shifts may take 
place. In moving upward along t^le diagram we find that the lowest 
layer in Figure 1 symbolically houses all personnel who are eithei: 
studying for the Bachelors Degree (subsector SI) or who have dropped 
out after some university or college training and are in industry 
(Subsector PI) . The non-degree holders layer allows for three types 
of personnel inflows. Thus, high school students can enter the 
student subsector via the HlSl path. Similarly, high school 
graduates from other nations can enter the same subsector via, the 
flow path labeled ElSl. The model allows for foreign nationals 
who have had some college or university training but who have hot 
attained the equivalent of a Bachelors Degree to enter the model 
via the path ElPl. It was a matter of convenience to allow these 
people to enter only into the professional (Pl) subsector. It is 
fully recognized that many of these people may, in fact, enter the 
country with the intention of studying further toward a Bachelors 
Degree and this may, in fact, be the predominant majority of such 
entrants. It was felt that in due time these people will, in fact, 
enter the student sector within the model. The availability of 
data, or better, the lack thereof, justifies this simplifying 
assumption in the construction of the model. Using similar 
reasoning, it was felt that all non-achievers who leave the country 
prior to obtaining their Bachelors Degree will leave via the 
professional subsector rather than leaving directly from the student 
subsector. The flow PlEl allows for all foreign nationals with some 
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undergraduiftte education to leave the country • The two remaining 
flows out of layer number 1 allow for the students who have, in 

fact, obtained the Bachelors Degree to move out and go on to anyone 

of the three subsectors in layer 2 which contains all personnel 
holding a Bachelors Degree. 

■Before moving into layer 2, however, we must account for all 
those people who leave layer 1 as the result of retirement or 
attrition. These are basically people who have completed their 
useful economic lives either as technicians or supporting personnel 
and are leaving the productive economy into permanent and complete 
retirement. The flow which allows for this to take place Is FlRl. 
Layer 1 allows for two intra-layer flows, specifically flow SlPl 

which channels all of the students who have dropped out into the 

professional or productive subsector PI. Similarly, flow PlSl 
allows for students who have at one time dropped out to reenter the 
academic or student subsector SI in order to continue their studies 
for the Bachelors Degree. It should be noticed that sub&ector PI 
represents that part of the economy, be it government, industry, 
commerce, or whatever, which employs people having some college 
training but no degree. 

The interlayer flow's can only come, in this model, from the 
student subsector of the layer .below. Thus, all those who have 
successfully completed their Baccalaureate studies leave subsector 
Si in layer 1 and move on to either of the three subsectors in 
layer 2. Layer 2, representing all those who do hold a Bachelors 
Degree but no more, differs from layer 1 in that withi.n this layer 
there are three subsectors. In addition to subsector S2, which 
represents all students studying for a Masters Degree, and subsector 
P2 which represents all those holding a Bachelors Degree but not more 
who are productively employed by our society in all areas other than 
university or college-level teaching, subsector F2 houses all 
Bachelors holders who hold no higher degree and who are employed in 
an academic capacity within higher education. 

Layer number 2 allows for only two inflows. Recent graduates 
with a Bachelors Degree from the layer below and foreign nationals 
.who enter the country holding a Bachelors Degree (E2P2) . Again, 
in this case as in layer 1, the assumption was made that all foreign 
nationals enter through the professional subsector » Similarly, all 
foreign nationals who leave this country holding a Bachelors Degree, 
but no more* leave through the professional subsector of layer 2 (P2) . 
In this case, retirement is allowed both from the professional 
subsector, namely P2R2 and also from the faculty or academic sub- 
sector, namely F2R2. It was felt that the retirment age of academics 
was generally somewhat different than the retirement age of people 
from industry or from government. Hence, the two flows were 
separated. 
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F i gure 1 

The Locations & Flow Paths of People with 
an Engineering Education: A Systems Schematic 

8 



In layer 2 there are not two but six possible intralayer 
flows. In addition to the population shifts between the students 
subsector S2 and the professional subsector P2, there are now 
similar paths for population shifts between either of these sub- 
sectors and the faculty subsector F2. All those who have 
successfully completed their studies toward a Masters Degree can 
leave layer 2 and move on to layer 3 via the channel labeled S23. 
Layer 3 is conceptually similar to layer .2 although Wa.thin it 
reside all those who hold a Masters Degree but no higher degrees. 
From layer 3 individuals may move up to layer 4 upon satisfactory 
completion of. the Doctorate Degree. Layer 4 is different from 
layers 2 or 3 in that within it reside all Doctorate holders and 
this layer does not contain a student subsector. Although it is 
recognized that there are post-doctoral students, it was felt that 
most of these people do hold- an academic rating, and so for the 
sake of simplicity, the student subsector was eliminated. 
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It should be noted that except for the exogenous inputs, all 
of the flows either originate and/or terminate at one of the sub- 
sectors. Therefore, &t each subsector it is, possible to write a 
"personnel" conservation equation in a manner analogous to the 
■coiiservation equations of either mass or energy or momentum in the 
physical sciences. These equations basically take into account, at 
any instant of .time, all, of the inflows, outflows and. rates of 
accumulation within a subsector. The equation states that the total 
rate bf flow into a subsector minus the rate of outflow is equal to 
the rate of accumulation .or growth. of population within a given 
subsector. Equation 1 states this in descriptive terms. It should 
be. noted that the rate of . generation term, quite often used in the 
physical sciences, is meaningless in .this case and therefore forced 
to equal zero. What we are basically saying here is that within a 
give?a subsector there can be no people generated. Any generation 
which does take, place in our population takes place some time prior 
to the individual’s entrance into either the academic or professional 
subsectors. Equation 1 is described in .mathematical symbols by 
equation 2. 

; Sum of , 

Input , 

Rates 



Sum of 
Output 
Rates 




u 



Rate • "~| 

Accumulation] 



( 1 ) 





( 2 ) 



Equation 2 may be rearranged from its differential formulation to an 
integral formulation which Is done in equation 3. 
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( 3 ) 



M 

t+dt t+dt 

/ [Jm - Jm Idt - / d(M) 

t 



If we now force the infinitesimal time increment dt to be equal 
to At a finite increment of time, and similarly, we allow the 
infinitesimal increment of accumulation dM to be equal to AM, 
then we may transform equation 3 to its difference format with but 
a few apologies to those who are inclined toward the purity of 
mathematics . 

Let 

dt -►At, and dM -► AM 

■ Vit • 

Allowing for some additional transformations such as those described 
by equation 5 

t 

If t - J, t -I- At « K, And At = DT ' (5) 

we can move on to equation 6 

- Mj + - iMjAt (6) 

and then on to equation 7 

M.K = M.J + (DT) - IMq] (7) 



which is in a format easily comprehended by the DYNAMO compiler 
developed in connection with a social-systems simulation methodology 
generally known as industrial dynamics which was developed by 
Professor J. W. Forrester at the Massachusetts Institute of Technology. 

Although the Idvel equations are determined by a straight 
forward application of the law of conservation of mass, the rate 
equations are determined largely by a group of relationships 
between many diverse factors. A brief summary of the "flowrates 
and the major influential factors" is given in Table 1. The table 
shows that the flowrates from one sector to another are determined 
by such gross influential factors as societal needs, group psychology. 
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student capabilities, national needs for faculty, student-faculty 
ratios, inputs, world technological gaps, and so forth. 

The approach that is used in combining these and determining t 
flow rate can perhaps be best explained by means of a number of 

specific examples. 



In Figure 2, we see an expanded version of the non-degree 
holders level where the student sector SI and the professional 
sector PI are shown in the square boxes. Let us consider the flow 
of professionals from Pi to SI and call that flowrate PlSl. In 
order to understand the dynamics of the flowrate PlSl we resort to 
a psychological submodel which is given in Figure 3. 

There are a numfcier of students studying for a Bachelors Degree 
in gubsector SI. After some length of time which we call a delay 
(say six months) we observe or count the number of students that 
exist in that sector. This delay is due to the fact that we do not 
usually count large number of students instantaneously, but it ta es 
some time before our bookkeeping system presents the data to decision 
makers . Furthermore , even after the figures regarding head count 
are available, it takes some additional time for' the decision makers 
to become really aware of what the figures mean. At the same time 
that there is an observation of the number of students, there exists 
in society a need for people to obtain a Bachelor s Degree namely, 

S need for students in the Bachelors program. This need is generated 
bv the discrepancy between the number of students who have a 
Bachelors Degree and those that we think we need with such a egree. 
It is reasonable to expect that if there is no difference between 
the need for students and the number of students that we observe 
there will be no dissatisfaction with the conditions that prevail. 
However, if we need more students than we have, there is s^e 
dissatisfaction generated. As the dissatisfaction grws there is 
a pressure to change the number of students that we observe, namely 
to open. the valve PlSl and let more students into school; 

The relationship between dissatisfaction and pressure to change 
the number of students is not a linear one. It can be seen from 
Fieure 3 that the relationship may be represented by an S-shaped 
curve. Thus, if there is very little dissatisfaction there is no 
pressure. But as the dissatisfaction grows, the pressure builds 
up rapidly. After the dissatisfaction has reached some large value, 
tL pressure does not continue to Increase linearly, but levels off 
and thereafter no matter how great the dissatisfaction, the pressure 
no longer increases. This kind of a relationship is representative 
of many real situations which we find in socio-economic situations. 
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Table 1 

FLOWRATES AND THE MAJOR INFLUENTIAL FACTORS 



Flowrate 

Professional ■> Student 
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Societal Needs 
Group Psychology 
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Student 
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Student Capability and/or 
School Standards 

Educational System 

(Student-Faculty Ratio) 
National Needs for Faculty 



Student ■> 



(student 
C Faculty 
(^Professional 



Structure of Educational 
System (Built in Delays) 
Current Level in Each 



s 

(Lower Level) (Higher Level) 
Foreign 

<j ^ Professional 

Personnel 



Higher Sector 

Past Inflows and Outflows 

of Lower Student Sector 

World Politics 

(Foreign Relations) 
Ecoilomics 

World Technological Gap 



Professional 1 
Faculty J 



Retirement 



Previous Graduation Rates 
Current Faculty and 
Professional Levels 



Need for 



Professional 



Doctorates 



Exogenous Input 

(Long Range Planning) 



High School 

> Student 



Historical Data 
(Extrapolated) 



Graduate 



A 










F i gure 2 

The Location i Flow Paths of 
Undergraduate Subsec tor ana 
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which "control" 
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F i gure 3 

A Typical Socio-economic Loop 
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In a similar manner , as the pressure Increases we tend to make 
some percentage change In the flow that we 'allow Into the student 
sector. Thlsi means that we might open the valve 2% or 3% greater 
than It was £it the previous Instant In time. The relationship 
between pressure and the percentage of change In the flowrate Is 
similar to the one given for dissatisfaction and pressure. From a 
mathematical standpoint we could combine the two functions Into one 
composite one, but this would obscure the physical significance of 
the variables and make it more difficult to make Intelligent 
modifications in the curve. We can see from Flsure 3 that the 
percentage change in flow to the student sector is the variable 
which will control the flowrate PlSl. Thus we have a dynamic loop 
in which the feedback, the delays, the need for students and some 
of the functional relationships that are postulated all conibine to 
regulate the valve. 

The actual variables used in this model are. shown in Figure 2. 

OSl represents the observation of the number of students that exist 
in the student sector SI at the present instant of time. This 
information together with the Societal need for students in SI, NSl, 
yield the variable DSl, which is the dissatisfaction with the number 
of students in SI. DSl provides information for generating pressure 
to change the number of students in SI, namely PCSl, and this in 
turn provides the information to generate CPIS, which is the 
percentage change in the flowrate. The change in the flowrate 
CPIS controls the valve PlSl. 

The mathematical relationships which represent the foregoing 
dynamic loop are represented by equations (8) through (15) . With 
the exception of a number of sumbols peculiar to the dynamo program 
and our own nomenclature, the equations are self~explanatory . It 
will be noted that- in equation 10, SID* (11. J) represents a special 
computer routine which stores the mmnber of students in SI at a 
given instant. in time and makes that number available sometime 
later (say 6 months later) in the analysis. In equations (13) and 
(14) , the symbols TPCSl and TCPS represent tables of values that , 
give the relationship between DSl and PCSl, as well as PCSl and 
CPIS. The tables are shown graphically beside the equations. 

Finally, the flow rate PlSl.KL which represents the flow rate during 
the time interval between time K and time L is given by the flow 
rate between time interval J and K plus some percentage of ? 
change in that flowrate. This. relationship is shown in equation (15). 



Sl.K - Sl.J + (DT)(H1S1.JK + ElSl.JK + PlSl.JK ; 
- SlPl.JK - S1P2.JK) 



08 ) 



Pl.K - Pl.J + (DTXElPl.JK + SlPl.JK - PlSi.JK 

- PlEl.JK - PlRl.JK) (9) 




PCSl.K - TPCSl (DSl.K) 



CPIS.K » TCP IS (PCSl.K) 



PCSl 



CPIS 




PlSl.KL - PlSl.JK + CPIS.K (PlSl.JK) (15) 



On the other hand, the flowrate from the student sector Si 
to the professional sector Pi is controlled by the valve SlPl. 

The flowrate is determined by some Internal as well as external 
factors. The external factors are similar to those given for the 
flowrate PlSl and are controlled by the external needs of our 
society for professional people NPl. The number of students who 
dropped out of school and do hot achieve a Bachelors Degree is 
determined not only by the need for professional personnel, but 
also by the mental ability of the students who enter the Bachelors 
pj^ogram and the educational standards of the school. The 
combination of student abilities and school standards cause some 
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percentage o£ the entering students to drop out. The percentage 
of students that drop out, PdIsI, can be assumed to be a constant 
or it may vary with time. For the purposes of the initial simulation 
we will assume that this dropout rate is a constant. The fraction 
of students ■ that will drop but due to the aforementioned causes, 
FRISI, is equal to the percentage of students who entered school 
from high school from foreign countries and from the professional 
sector. The studehts do not drop but all at one time. Mostbf 
them drop out during the first yeat bf cbllege and less and less 
students drop out as the years go by. This phencraenon is taken 
Into accbuiit in the. model by introducing a special routine known 
as a delay. The delay routine takes the flowrate and distributes 
it in an appropriate manner over the time period involved. The' 
flowrate SlPl is thus the sum of the internal flowrate ESIPI. Some 
of the definitions of terms arid typical equatibiis that are used in 
the liitemal loop are given below. 

The rate at which students graduate from SI and who proceed 
to level 2 is called S12. This flowrate is determined by the rate 
of students entering as freshman * approximately 4 years earlier than 
the present time minus the rate at which students dropped out over 
the 4 year period* All bf these rates are delayed in an appropriate 
manner by a delay routine. 



er|c 



PDISI = Constant or Variable with Time 

* Percent Dropout of Students Entering SI 

Diie to Mental Abilities, i.e.. Internal Reasons 



FRISI.KL = PDISI (HlSl.JK + ElSl.JK + PlSl.JK) (16) 

ISIPI.KL = Delay 1 (FRISI.KL, DRDSl) (17) 

DRDSl * Cbnstant 

“s Dropout Delay at SI in Time Units 
SlPl.KL = ISIPI.KL = ESIPI.KI, (18) 

The rate at whibh people retire from the prbfessional sector 
PI is called PlRl. As might be expected, that flowrate is determined 
by the number of people that entered this system 30 or 40 years 
earlier from high school, the number of pebple who retire from the 
Masters and Ph.D. sectors and some , appropriate distribution functions 
regarding the age at which people retire. 



Finally, the rate of inflow and outflow of foreigners from 
the professional sectpr, namely ElPl and Fl£l are recognized 
explicitly in the model and their respective values are assumed to 
be reasonably constant. An extensive search of available 
information regarding foreign personnel entering and leaving the 
United States has revealed that very little good data is now being 
kept on this subject. If we were to assume relationships which 
are other than constat values for these flowrates , we would have 
to greatly expand the model to Include the Influence of foreign 
economics, politics, world affairs, cultural patterns, and so forth. 
For the present time we will make reasonable assumptions regarding 
these flowrates and turn our attention to more realistic 
assumptions at a future date. 

In a manner substantially similar to the one presented for 
level 1, the model has been developed for the second, third and 
fourth level. The flow diagrams are shown in Figures 4, 5 and 6. 



SUGGESTED FUTURE EXTENSIONS OF THE MODEL 

At the present time the model and the simulation are concerned 
strictly with the profession of Engineering (not distinguishing 
between mechanical, electrical, etc. engineers) . However, in the 
future we are looking forward to doing the following. 

a. ND4 in the model represents the need for doctorates 
which is felt in the general population at a given 
time. This variable mus.t be functionally related to 
the projected long-range needs of the nation.. An 
extension of the present model is required. 

b. Applications to other professions and disciplines each 
treated individually. 

c. Generalization of the model to in 9 lude several inter- 
related professions and/or disciplines. (Physical Sciences, 
Mathematics, and Engineering, The Life Sciences , Behavioral 
Sciences, and Medlcinej The Social Sciences, Law, and 
Business Administration . ) 

d. Extension of the generalized model to interrelate all 
areas of higher learning. 

e. To further delineate within, all of foregoing model the 
impact of economic variables (i.e. , scholarships., 
industrial, government and academic salaries, avaiability 
of money for educational plant...) 
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F I gure 4 

. . . c I nuu Put-h<^ of People within the 

^Polt'eacJaliureate ^iJbsector and the 

ihich '’control" the Population Flows 
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FiCjUre 5 

The Location & Flow Paths of People within the 
Post Masters Subsector and the Interrelation 
of Socio-economic Factors which ’’control’' 
the Populat i on Flows 
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Figure 6 



The Location i Flow Paths of People within the 
p5It Doctorate Subsector and the Interrelation 
of Socio-economic Factors which control 



the Populat i on Flows 
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APPENDIX 



NOTE p rint cards START AT AOOO O 



NOTE 

SPE^ 

NOTE 

2L 



"DEC. 27 » 1967 
nT=l /I FNCTH=010/PRTPER=2/PLTPER=2 



.S K K= SI. J + ( DT)(HlSl.JK-t-ElSl.JK-t-PlSl.JK-SlPl.JK>Sl2.JK:^ 



7N 

2L 



8N 



pl,K=Pj ._J -^nTi ( ElPl.JK-HSlPl.JK-PlEl. 

"pT^r2‘+P3+‘PA 

H1S1.KL=H1S1.JK • ^ — __ ^ 



6N 



H1S1=62000/12 

E1S1.KL=E1S1.JK 



E1S1=4600/12 

1 4R P ISl . KL=P1S1. JK-h(C'P1 S.K) ( PI JI^JJO. 

“6N PlSi=0 

3-7B* S1D=B0XLIN(7, 1) 



PI "TO SI 



OSl.K=S10>!'7.K 
S1D*1.K=S1.K 



SlD=BOXLOAb(Sl, 1) 

6N 0S1=S1 J. 

T6"A *NST.K=NS2.K+NF2.K+NP2.K-F2.K-P2.K-S2.K 

6N NS1=0S1 

DST7i^NSl.k70Sl.K 

6N DS1=1 



NE"ED F.OR SI 



P(:S1.K=CLIP(PCS1P0S.K,PCS1NEG.K,DS1.K|1) 

PC sipos. K= 1 -EXP ( (-A) ( ( D^i. K-1 ) ~ ^ ^ jj — — 

ptTTw F'n7K ^^^pT(“4TT^^ . k1 TDS 1 .K)(DS1.K)) 

CP1S.K=(PCS1.K/P CSIAB S.K) {1-EXP( (-A)(P C S1ABS.K)-PC SU_ BS.K-PCS1^ 

Tk“pc sTaFsTk") ) " ■ 

PCSIABS. K= MAX (PCS1.K,-PCS1.K) 



IN act. matr. for bach, student^ 



TH TTZ= 20200 / 1 2 

DMS1=45 DEL_A^ . 

"18R — TiS 1 E“XJ^'( S"UMS . JK ) ( 1-ADOS 1 . K ) 

SD0S1.K=SUM1( 52,RAVS1) 

DROP OUT OF SI 



44A 



ADOSl.K=SbbTl. K/52 AV 
BAVSl^l. K= ( S1 P1.JK)(DT)/S1.K 
"BAV S'l = BOX L TN ( 5 2 , 1 ) 

J U M S . ia^BSMSl »52 .JK 



“8 R B S fTS~l * 1 . k"L= HIS 1 . J k + E1S1.JK+P1S1.JK 

B SM S l = B OXLINj 

BSMSl~=BbXLbAD( HlSl+ElSl+Pl SI 1 1 ) 

36N RAVS1=B0XL0AD( BAVS1»1, 1 ) 



7A 



DTPY=12 MONTHS' 
S1P1.KL=IS1P1. JK-hESIPI. JK 



6N 



MOTE 

6N 



S1P1=8200 

IS1P1.KL=DE LAY3(FRIS1.JK> 2 A) 

PRECEDTNG S“CH0ULD be DE'LAY 1, NOT 3 
IS1P1=0 



T9”R FRIS1.KX= (PDISl) (HlSl. JK+EISI. JK+PISI. JK+0 ) 

PDIS1 = 20 9/12 ( . 5-Sl ) /1 2 

TaF ES1P1.KL=ES1P1.JK+( CS1P.K)(ES1P1.JK) 

6N ES1P1=0 



SOCIAL EXTERNAL 



PCP1.K=CLIP(P"CP1P0S.K,P-CP1NEG.K,DP1.K,1) 

pcpiPOS.K=i-Exp( (-4) ( (DPi.K-1) L^P ij Jiri-U ^ ^ * *Szi-LLL 

— — =rrr-rr — rwT-i i/_nni i/ \ \ 



PCP1NEG.K=-EX”P ( ( -A) ( DPI. K-DPl. K-DPl. K) ) 
PCP1ABS.K= MAX (PC P 1_._K -PCP I.K) 



‘CSTp7k= (PCPI.K / PCi'P l^ABS . K )'( 1- EX P ( ( -A ) - PC PI PWE R .K ) ) 
PCPIPW ER. K=PCP1A 
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6N 

‘7A‘ 

7N 

"12A 

7A 



■ O pT; K^NPT rKyO PI . K“ 
DP1=0 



28A 

29A 

c 

C 

C 

12N 



’’NPI" k=NN 1 . k-N SI . K 

NP1=NN1-NS_1 __ 

wr.'k="l DMi X f . K ) ( ND4. K ) 

dm I X 1 . K= CC l_t.AT_B 1 . K 

XfBlTK^'l AA’l ) E"XP I R LT 1 . K ) 
BLT 1 . K = ( BB 1 )LO GNMTIM E.K) 
AA1=0 

BB1=0 



needs OF PI 



6N 

C 



'CC1=17 

N N 1 ^ D MIXl) (ND^) 
"DMIX1=20 

N D4=650/ 12 

0P1.K=P1D*7.K 
PlD’i‘1. K=P1.]^ 



■pfD=BOXLIN(7, 1) 
PlD=BOXLOAD(Pl t 1) 



bR _ 
6N 
9R_ 
■6R“ 
19N 



0P1=P1 
PlRi,K L=RRl » JK 
■pfRl=0 



RRl.KL = TOTR ^-RR 2 . J K-RR 3 . J K-RR J JL. 



6R 
6N . 
6R 
6N 
2L 

_2_^ 
8A * 
2L 
20N' 
_2L_ 
2 ON 

JA 

44R 
44 R 



foreign to pi 



PI TO El 



8R 
20N 
T4R ■ 
6N 



"T0TR.KL= TOTR. JK 
TOTR^ISJ (^1+P 2+P3_4 -P4j. 

TrPl.KL = ElPl.JK 
E_l£l^ - 

'PlEloKL=PiEl. JK 

-r^^!2TmmiTi2;TK+T2S2nK“FTs"2^^ ■ 

S2= 12400 /DTP Y — nE""p7 

“Asis2?i=r2lK-Fl^.K 

lTS27kL=IS12.JK) (S2.K)/AS1S2.K 

S 1 P 2 . K L = I S 1 2 » J K ) ( P 2 . K ) / 

"SlF2.kT=S12.'jK-SlS2. JK-S1P2.JK 
S1F2=1000 /DTPY 



F 2 STTkr^lTfrjK + (CF2S.K)(F2S2.JK) 
F2S2=0 



PCS2.K=TLIPlPCS2PdrrK,PCS2NEG.K,bS2.K,i) , . 

PCS2P0S.K=1-EX P( ( -^ ) ( ( D S 2 » K- 1 M P_SA*JilliJL9. ^ ^ — 

prs'2NEG'. K^^^EXTI ( - 4 ) { DS 2 . K- D S 2 . K- DS 2 . K ) ) 

^• 72 l 7 K 7 ;K^-^:{ypcl 2 A^ 

PCS2PWER.K= P C_S_^AJ S i A BSj^^ 



TOA FjT2~.'I^N^» K+NF3. K+NP3. K-F3."K-P3. K-S3.K 

6N NS2=0S2 

D'S2.K=MS2.k/0S2.K 

6 A 0S2.K=S2D’!'7^IA 



NEEDS OF S2 



6 MONTH DELAY 



“>2D’.n.K=S2.K 
S2D=B0XLIN(7 t 1) 



6N 



S2D=BdXL0AD(S2, 1) 
0S2=S2 
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lAR 

20N 



8A 



29A 

8N 

6N 

C 



C 

C 



6A 



lAR 



P 2 S 2 «KL-P2S2»J K + (CR2S«K) (P2S2« JK) — 

P2S2=0 

CP2S.K=(PCS2.K/PC S2ABS.K) ( 1- EX P( (-4 ) ( PCS2PWER «K) )1 
TzTz. kI= S2 P 2 . J K + (, CS2 P K ) ( S2 P2“.J K ) 

S2P2=6500/DTPY_ 



~CS2’P.K=(’PCP2rK/PCP2ABS.K)(l-EXP(-4-PCP2PWER.K).) 

CF2P.K=CS2P«K 



lAR "F2P2.'KL=F2P2. JK + (CF2P.K) ( F2P2.JK) 
6N F2P2=0 



8A 

20A 



PCP2.K=CLIP(PCP2P0S.KtPCP2NEG.K,DP2.Ktl) ^ ^ 
PCP2P0S.K=1-E XP( (-A^) ( (DP2.K-1) (DP2.K-1) (DP2«K-l )J^j_ 



PCP2NEG. K=-EXP(-A-( DP2. K-DP2. K-DP2. K) ) 
PCP2ABS.K=MAX( PCP2.jViP^P2jJ^ ' 



"P'C'P~2 PWi R • K=P CP 2 A'BS . K-PC P2 A BS . K- PC P2A BS . K 
NP2.K=NB2«K-NF2«K-NS2»K 



6N NP2=0P2 

12A NB2.K=(DMI X2.K ) (NDA,K) 



TA DMIX2.K=CC2+ATB.K 
2 8A ATB . K= ( A A 2 )E XP(BLT.K) 



"B IT . K= ( BB2 )T0 GN ( f I ME . K ) 
NB2=NP2+NF2+NS2 



DMIX2=76.3/5.00 
AA2=0 



BB2=0 

CC2=13 



DP2.K=NP2. K/0P2.K 
0P2.K=B0XP2*13.K 



6A • B0XP2*1.K=P2.K 
37B B0XP2=B 0XLIN ( 1 3t 1) 



“B0XP2=B0XL0AD(P2t 1) 
cyF7.KI=<;?F2«JK + (CS2F.K) (S2F2« JK) 



«2F^K= ( PCF 2 .K/PCF2AB S.K) ( 1-EXP (- A-PCF^WER_^) 
PCF2PWER. K=PCF2ABS. K-PCF-2ABS. K-PCF2ABS .K 

PCF2ABS.K=MA X ( P CF2 . Kt -PCF2 .K ) 

~7A PCF2.K=PCF’2E. K+PSFR2.K 

7N PCF2=PC F2E+PSF R2 



p's FT2 Tk^CT IP(PSFR 2P0S.KtPSFR2MEG.KtSFR2.KtMSFR2.K) 

PSFR2P0S. K=l-E XP (-A-SFR2P0V<ERABS« K ) 



SFR2P1^'R» K= ( SFR2.K-MSFR2) ( SFR2.K-MSFR2.K) (SFR2.K-MSFR2) 

S F R 2 P0WE RABS.K=MAX (SFR2PPER.K t-SFR2P^ .J<J ^ 



■pSF rTNEG . K=- e XPT- A-S FR 2 . K-S FR 2 . K- S F R2 . K ) 

MSFR2=20 AVE R. OR DESIRED MASTERS LEVEL START OF SFR2_ 

Tf.K=F2.K+F3. K+FA.K . 

SFR2.K=S2«K/FF,K 



8N FF=F2+F3+FA 
20N SFR2=S2/FF 



PCF2E .K=CLI P ( PCF2EPOS. Kt PCF2EMEG. Kt DF2.K»1 ) 

PC F 2 E PO S • K = 1-EXP ( -A- ((DF2«K-1) (DF2 «K-1) (DF2«K-1) )) 



19A 

6N 

7A 



PCF2ENEG.K=-EXP f-A-C DF2 . K-DF2. K-DF2. K) ) 

DF2.K =NF2«K-nF2« K 

"dF2=1 

NF2.K=(FMIX2.-K1 ( NFA. K+NF3* K +NF2* KjtPJ- 
NF2=OF2 

FMIX2.K=FF2+DT E«K 



NEEDS OF F2 



28A DTE.K=(DD2)EXP(ELT.K) 
29A ELT . Kf±EE 2 ) LOGN (TIME . K ). 
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ERIC 

^skisqeb 



c_ J E_2=p 

C FF2=.l ■ 

j, A OF 2 . K= B01FJ«l.j<: 

6A BOXF2«i.K=:F2.K 

B 0 XF.2 =BOX UNI?. iJJ— 
B0XF2=BbxL0AD( F2t 1) 
6N OF2 =F2 



]^_NIH_DELAY. 




lAR 

6N 



5 OR 



P2F2iKL=P2F2.JK+(CP2F.K) ( P2F2.JK) ^ 

£2j^2.ja5JAfL2.JJi^^^ 



20N 



P2R2=2AAA/bfPY 

R R 2 . K L = P E la y 6 j J XS ^*5 LLtlftl 
”^Si2=BbXLfNl516t 1 ) 

B X S 1 2^ boxlo ad ( S_l_2 i 1 ) 



6A "bXS12«"1.K=S12.v 

P7D7.KI = RR2. JK-P2R2* _JJ< 
T2R2=156/DTPY ‘ 

S 2 3 • K L = DE LAY 3 ( MS2EX* JK iAZ 

’"shoulIT have been delay 1» 

S 23=AA AA/DTPY 

TBR “MS2EX.KL= (SUMA. JK) ( 1-ADOS2.K) 
S UMA.K= BSMJ.2*Jl5jl]1 



20N 

"notF 

2 ON 




8R 

37B 



S UM A » N= p c6 '~lyt 

BSMS2’!‘l. KL=F2S2. JK+P2S2. JK+S1S2. J 

B S M S2=BOX LIjiLF5jLll 



53 A 

22A. 

‘37B 



BAVS2=BOXLIN(35» 1) 

B S M S2=BOXLQAD ( F2S2+P2S j. +SlS2) 1 J_ 
’^VS2=BOXLOAD( BAVS2*lt D 



6R E 2 P 2 » K L = E 2 P 2. »_JJi. 



20N E2P*2=3047/DTPY 

P2E 2 .KL=P2E2, JK_ 



P2E2=0 



note layer 3 

-2L ^Tn^:j^(0T)(S23.jKtP3»~.JK^F3S3.JK~rs34.JK-S3P3.JK- S3F3.JK) 

2 ON S3 = 292 1/LJPY 



ir rS2ljl^]:K_.P3,K.F3,^^ 



8A 

2L_ 

20N 



A4R 

A4R 



8R 

lAR 



l#3\"Kr=Tr2F: J K )T S"3^^M S3 . K 

PC S3 NEG . K= -E X Pjr Jtri D-il 



IS S2S3 CORRECT 



_^_S3 A BliJ!L= M AXJ P CL3 ._K»jiPCS3 •.!<-» 
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• 


DS3.K=NS3. K/0S3.K 




8A 

6N 


NS 3 . K= ND4 . k- P4 . K-F 4 . K 

NS3=0S3 




'6 A 


DS3.K=SlO’^7.K 6 MO DELAY 

■S3D=«‘1 *K=S1.K 




37B 


S3'0=BbXLIN(7tl ) 

S^nsROXI HAni S3. 1 ) 


' 


-NOTE 


bS3=S3 

DELS3=48 WHATS IT FOR ^ 




-lAR 

6N 


"P3S3.KL=P3S3.JK+ICP3S.K) (P3S3.JK) 

■p3S3=0 




14R 


CF3S-.K-CP3S.K SAME FOR NOW 

S3P3.KL=S3P3.JK+|CS3P.K)(S3P3.JK) 




6N- . 


S3P3“0 

rS3P - K= ( PCP3.K/PCP3ABS.K) ( 1-EXP( -4-PCP3PWER .K) ) 


’ 


PCP3.K=CLiP(PCP3P0S.K,PCP3NEG.K,DP3.K,l) 

•prP^pn<;.K=1-EXP(-4-( (DP3.K-1) (DP3.K-inDP3»K~_l_nj 


“ 


PCP3NEG.K=-E XP (-4-( DP3.K-DP3.K-DP3#K) ) 

• prpapv.|FR.K = prP3ARS-K-PCP3ABS.K-PCP3ABS.K 




PCP3ABS.K=MAX( PCP3.Kt-PCP3.K) 

bP3;K=NP3.K/nP3.K-- ■ ■■ — 


0P3.K=B0XP3*7.K' 

BOXP3*l.K=P3.K 


37 B 
6N 


B0XP3=B0XLIN(7tl) 

0P3=P3‘ 




36N 

8A 


BbX'P3=B0XLbADI P3tl) ■ 

NP3.K=NM3.K-NF3.K-NS3.K 




6N 

12A 


NP3=0P3 

NM3.K=IDMIX3.K) (ND4.K) 


• 


7A 
2 8A 


DMIX3.k="CC3+ATB3.K 

ATB3.K=I AA3)EXPIBLT3.K) 




29A 

8N 


"”B LT 3 . K= 1 B B 3 ) LO GN I T I M E . K ) 

NM3=NP3+NF3+NS3 




6N 

c 


DM 1X3= 18. 8 15 

AA3=1 • ' • 




C 

r 


BB3=1 

r.r.3=?o. 0/6.0 




14R 
. 6N 


F3P3.kL=F3P3,JK+(CF3P.K)(F3P3.JK) 

F3P3=0 ^ 


• 


14R 


CF3P.K=CS3P.K SAME FOR NOW 

S3F3.KL=S3F3.JK + ICS-3F.K) IS3F3.JK) ^ 




6N 


S3F3=0 

CS3F.K= IPCF3. K/PCF3ABS.K) 1 1-EXP 1 -4-PCF3PWER .K)_). 




7A 

7N 


PCF3.K=PCF3E.K+PSFR3.K 

PCF3=PCF3E+PSFR3 




m 


PCF3ABS.K=MAXI PCF3. Kt-PCF3.K) 
PCF3PWER.K=PCF3ABS.K-PCF3ABS.K-PCF3ABS .K 






PSFR3.K=CLIPIPSFR3P0S.K,PSFR3NEG.K,SFR3.Ktl) 
PSFR3P0S.K=l-EXPI-^4-(SFR3.K-l) ISFR3.K-1) ISFR3.K-1)) 




26A 


PSFR3NEG.K=-EXP|-4-SFR3. K-SFR3. K-SFR3 . K ) 
SFR3. K= ( S3'. K+0+0 ) /( F2. K+F3. K+F4.K) 




. 26N 


SFR3=(S3+0+0)/IF2+F3+F4) 

np:».l^=MF3.K/nF3.K 




PCF3E.K=CLi P (PCF3EP0S.Kt PCF3ENEG. Kt DF3. Kt 1 ) 

PCF3EP0S.K=l-EXPl-4-inF3.K-l) IDFS.K-niDFLiKzIJi 


• 


PCF3ENEG.K=-EX’PI-4-0F3.K-DF3.K-DF3.K) 
NF3.K=IFMIX3.K) (NF4.K+NF3.K+NF2.K) 





er|c 



28A 

‘29A 

6N_ 

.... 

C_ 

"C 

6A 



"FMI X3 . K=Tf T+DTE 3 . K 
DT E3 . K= I DD3 [EX P I.ELJ.?* J<1- 
■ ELT3 . K= i EE3TLDGM I TI ME. K ) 

NF3=D F3 

■'dD3=T 

EE3=0 



6A 

J7B_ 

6N 

36N 



"FF3=0.3 

DF3.K= B0XF3« 7.K 

■BDXF3«1.K=F3. K 
B D X ^ = B D XLIN(7tl) 

D r3=FT' 7 

B DXF3= 3 D X LD a PlFJj.lj] 



39R 
_MDIE_ 
Ton 
^R_ 
'5 3 A 
_20^ 
22A 
37B 



r3 TTK L =d“E L AY 3 I MS 3 E X . J K , 1 2 ) 

S HDU UL B E D ELAY It NOT 3 

S3A=600/DT"PY ■■ 

MS 3 E X . K L = ( S UMB . JK)(1-ADDS 
’/ruyrD.T= S0MT62 1 BAVS3. K ). 
ADDS3.K=AUX A P. 



6A „ 
"7 A 
37B 



AUUO^« I , .1 1 f I-IT \ \ 

Wsl#iTK=Tl/S3. K) I I S3P3. K) I DT ) +1 S3F3. K) IDT)) 

BAVS3=BDXLIN I62jj_) : 

BTTv'SBTBOXLDAb I BAVSB*! i 1 ) 

SlJMB.K=BSMS3Ji52J< ; 

TSMS"3"«r. K=P3S3. JK+S2S3. JK 
_BSJiS3 = BD X_L I M.1.5 2jJL) 



50R 

39R_ 

'39R 



BS7^S’3=BWLDAb|P3S3+S2S3t 1) 

P 3 R KL=(RR3.JK)(P3.K)/( P 3 . K Kj.. 
p"3R3=‘ 



RR3 .KL=DEUW3I RR3A._JKj.p_R.EI3L.. 
■“TR3A'.TL = bE LAY 3 1 UDRR3 . J K , DRET3 ) 



44 N 
C 

AA 


nRET3=(UYK3J lUIKT 
DYR3=l2.b) IDtPY) 
iinppr; ,k = RRR3*31 . K 




OM 

6A 

^ “7 D 


B'RR3*1.K=S237j'k' 
noDR-Rfiyl TN(31.SFTR3) 




3 > D 

A.4N 


“B RR3 = BD X LO A b I S2 3 1 1 ) 
QPTRRs ( YR3) 1 DTPY) /30 


— 


7N 

r 


“YR3=YR3T-2 ^ ^ 

YR3T-3B TDTAL YEARS 


TO RPTTRFMENT AFTER MASTER^PIP 


w 

• 7R 
1 ?N 


■ F3R3.KL=RR3.JK+P3R3.JK 
E3P3-I.1) IS23) 


1 + DR - ) 


6N' 

AN 


E3P3.KL=E3P3. JK 
DRFR.KI =P3E3.JK 


f 


OIN 

1 4R 


P3E3=b 

P3PR,kI=PRF3-JK+ICP3F.K)|P3F3.JK) 




7N 


"P3Fr=PC“F'3E+PSFR3 

CP3F.K=IPCF3.K/PCF3ABS.K)I I-EXPL-A; 


-PCF3PWER.K )) 



NDTE 

2L 

12N 



C__ 

”2L 

20N 



50 R 



P(^p«7omtTrn?CTirrFTf^^ .jk-p4e<. . jk-p4r«. jk) 

»=I PC TP A) (D4I ) — 

VI=9000/DTPY°°^* "^P.P.CTPLATE. 

4Vk = F 4'. j +T DT M S ‘3 F 4 . Jk + P 4 F 4*. J K- F 4 P'4 . J K - F 4 R 4 . J K -tO +0 ) 

!^=2720/pTPY — 

3F4.K=S34. JK-S3P4.JK 

3P4.KL=IS34.JK)|P4 ..K.)_/ ,I_P 4 




28 



14R F4P4.KL=F4P4. JK + (CF4P.K) ( F4P4. JK) 

% # n — HTTlFr#"/.— nr da DU CD _ 



er|c 



TF4P7k=7pCPA. k7 PCPAABS. K) ( 1-EXP"( 4-pTP4PWER .K ) ) 
0p/t,K=NP4,K/nP4.K 



P C74Tk= C U P ( p CP.4 Pd S. K , P C P4N E G'. K , D p7. K , 1 ) 

PCP4P0S.K=1-EXP( -4 -(DPA.K -l) (DP4.K-1) (DP4.K-1 )) 



P"CP4NEG.K=-EXP(-4-DP4. K-DP4, k-bP4,K) 

PC P4PWE R.K=PCP4A BS.K-PCP4ABS. K-PCP4ABS. K 

TCPAABS. K=MAX rPCP4,K,-PCP4. K) 





6A 

6A 


UrH • DU A r 1 j • IN — 

bOx”Pa*i.k=pa. k . . ■ ^ 

BOXPA = BOXLIN( 13f 1) ^ 


B0XP4=B0XLC)A0(P4, 1) 

nF4.K=NF4.K/0F4.K . 


n;|-p4.k=Rr)4.K-NFA.K 

ND4.K=D4.K ^ — : 


U4.K=U4.J 

D4=F4 = 




NOTE 

6A 


04 VALUE UbtO y NOT CORRECT* 

0FA,K=B0XFA*7. K 6 MO DELAY 




6 A 
37B 


BdXFA«i.K=FA. K 

B0XFA=B0XLIN (7t 1) 




36N 

6N 


BdXFA=BOXL0AD(FA, 1) 

np4-P4 




f9A 

^N 


"NFA.K=TFMiXA. K) ( NFA. K+NF3.. K+NF2. K+0 ) 

MP4=nF4 — 




8A 

lAR 


FM I XA . k= i-FM I X3 . K-FM I X2 . K 

PAFA.KL=PAFA.JK+(CPAF.K) (PAFA. JK) 




6N 


PAFA=0 

pr,F4.K = ri TPf PCFAPOS.KtPCFANEG.KjDhA.K, 1)_ _ 


' ’ PQP4PQS. K=1“EXP ( “A-( DF4. K-1 ) ( DF4. K-1 ) ( DFh.K 1)) 

PCF4NEG.K=-EXP(-4-0F4.K-0F4.K-DF_4.K) ___ 


■ " " PCT4aB"S.K=KAX(PCF4.K,-PCF4.K) . ^ 

PCF4PWER.K=PCF4ABS .K-PCF4ABS. K-PCF4ABS .K 


roAC i^-/DrF4.K/PCF4ABS.K)‘( l-EXP(-4-PCF4PWER.K )) 




NOTE 

NOTE 


■p'ARA. kL=(RRA. JK) rPA.K)/(PA.K+FA.K) 

P4R4 CHANGED BECAUSE OF SIMUTANEOUS EO. 




ION 


P4RA.KL=RRA.JK-PA.k 

PARA=S3PA+EAPA-PAEA-PARA+FARA+0 




. 7R 


FARA. KL=RRA. JK-PARA. JK 

F4R4=0 - — 




NOTE 


ABOVE IS INCORRECTf CHANGED TO RIO SIMULTAN EO. 
RR4.KL=DELAY6(BRRA=:=25,2A) 




6A 

37B 


B'RRA*1.K=S3A, JK ■ 

BRR4=B0XL IN ( 21 » 1 ) 




12N 


BRRA=BOXLOAO( S3At 1) 

RRA=(. 15M9000) 




6R 

^N 


PAEA.KL=PAEA. JK 

PAEA = 0 




6R 

6N 


EAPA.kL=EAPA. JK 

E4P4 = 0 




NOTE 

KiriTP 


ENO OF OYNAMO DECK 

PRINT CARDS . — . — TT — L-TTrr-r;r::rV 


PRINT ' l)Sl»S2f S3/2)P1»P2»P3» PA/3) F2» F3, FA/A) SI 2 » SI S2 t S 1P2 1 SIPZ/ P » i 

,„lt,;,:?"F,/MS3«.»,S3P4,S3F4/7)SlPl,S2P2,S3P3,»/S)PlSl,P2S2iP3S 




X 

PLOT 


^*S1 = S.P1=P/S2P2=B, F2=C/S3=D, P3=E_t.F3=F/ P^'=Gf.FA=H 
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